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Starting from the Liouville equation, and using a BBGKY-like hierarciiy, we derive a kinetic 
equation for the point vortex gas in two-dimensional (2D) hydrodynamics, taking two-body corre- 
lations and collective effects into account. This equation is valid at the order 1/N where N ^ 1 \s 
the number of point vortices in the system (we assume that their individual circulation scales like 
7 ^ i/N). It gives the first correction, due to graininess and correlation effects, to the 2D Euler 
equation that is obtained for A'^ +00. For axisymmetric distributions, this kinetic equation does 
not relax towards the Boltzmann distribution of statistical equilibrium. This implies either that (i) 
the "coUisional" (correlational) relaxation time is larger than Nto, where in is the dynamical time, 
so that three-body, four-body... correlations must be taken into account in the kinetic theory, or 
(ii) that the point vortex gas is non-ergodic (or does not mix well) and will never attain statistical 
equilibrium. Non-axisymmetric distributions may relax towards the Boltzmann distribution on a 
timescale of the order Nto due to the existence of additional resonances, but this is hard to prove 
from the kinetic theory. On the other hand, 2D Euler unstable vortex distributions can experience 
a process of "collisionless" (correlationless) violent relaxation towards a non-Boltzmannian quasis- 
tationary state (QSS) on a very short timescale of the order of a few dynamical times. This QSS 
is possibly described by the Miller-Robert-Sommeria (MRS) statistical theory which is the coun- 
terpart, in the context of two-dimensional hydrodynamics, of the Lynden-Bell statistical theory of 
violent relaxation in stellar dynamics. 



I. INTRODUCTION 

In 1949, Onsager [l[ published a seminal paper in which he laid dovifn the foundations of the statistical mechanics 
of vortices in two-dimensional hydrodynamics [s^ . He considered the point vortex gas as an idealization of more 
realistic vorticity fields and discovered that negative temperature states are possible for this system [8^1 • At negative 
temperatures, corresponding to high energies, like-sign vortices have the tendency to cluster into "supervortices" 
similar to the large-scale vortices (e.g. Jupiter's great red spot) observed in the atmosphere of giant planets. If the 
point vortices all have the same sign, one gets a monopole. If they have different signs, one gets a dipole made of two 
clusters of opposite sign, or possibly a tripole made of a central vortex of a given sign surrounded by two vortices of 
opposite sign. Therefore, statistical mechanics may explain the ubiquity of large-scale vortices observed in geophysical 
and astrophysical flows. 

The qualitative arguments of Onsager were developed more quantitatively in a mean field approximation by Joyce 
& Montgomery 0, Q] , Kida Q and Pointin & Lundgren Q , and by Onsager himself in unpublished notes [1] . 
The statistical theory predicts that the point vortex gas should relax towards an equilibrium state described by the 
Boltzmann distribution. Specifically, the equilibrium stream function is solution of a Boltzmann-Poisson equation. 
Several mathematical works have shown how a proper thermodynamic limit could be rigorously defined for the 

point vortex gas (in the Onsager picture). It is shown that the mean field approximation becomes exact in the limit 
N — !> +CXD with 7 ^ 1/N (where N is the number of point vortices and 7 is the circulation of a point vortex). 

A practical limitation of Onsager's theory resides in the approximation of a continuous fiow by a discrete collection 
of point vortices. This is clearly an idealization that he was aware of: "The distributions of vorticity which occur 
in the actual flow of normal liquids are continuous" [l|. This approximation also leads to some ambiguity in the 
construction of a statistical model of realistic flows: "What set of discrete vortices will best approximate a continuous 
distribution of vorticity?" Q . 

The statistical mechanics of continuous vorticity fields was developed later by Miller [T^l and Robert & Sommeria 
[l5| . The same theory appears in a paper (in russian) by Kuzmin [l6| eight years before but his contribution is less 
well-known. The Miller-Robert-Sommeria (MRS) statistical theory is based on the 2D Euler equation which describes 
incompressible and inviscid flows. It predicts that the 2D Euler equation can reach a statistical equilibrium state (or 
metaequilibrium state) on a coarse-grained scale as a result of a mixing process. Recently, the MRS theory has been 
applied to _geophysical and astro phy sical flows, notably to oceanic circulation [l3, jovian vortices (Jupiter's great red 
spot) [H-io], Fofonoff flows [2ll42a| and oceanic rings and jets 

The MRS theory shares many analogies with the theory of violent relaxation developed by Lynden-Bell [2^ for 
collisionless stellar systems described by the Vlasov equation. The analogy between two-dimensional vortices and 
self-gravitating systems was mentioned by Onsager in a letter to Lin: "At negative temperatures, the appropriate 
statistical methods have analogues not in the theory of electrolytes, but in the statistics of stars..." 0. This analogy 
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has been further developed in [261 ]. 

At that stage, we would like to make some observations that will be important in the following. First of all, 
point vortices do not exist in nature [IBl- Fundamentally, the physical system to consider is the continuous vorticity 
field that is solution of the 2D Euler equation. The point vortex model is a particular solution of the 2D Euler 



equation of the form uj{r,t) = J2iLiliH^ ~ i'i(^))- However, there are rigorous results [27[ which show that any 
smooth solution a;(r, t) of the 2D Euler equation may be approximated arbitrarily well over a finite time interval by 
a collection of point vortices with 7^ ~ 1/N as N ^ +00. The physical situation is completely different in the case 
of stellar systems. Basically, a stellar system is a discrete collection of N point mass stars. For N +00 in a proper 
thermodynamic limit, this discrete system can be approximated by a continuous distribution function that is solution 
of the Vlasov equation. This is the reverse situation with respect to 2D hydrodynamics! It is nevertheless interesting 
to further develop the analogy between a system of N stars in astrophysics and a system of N point vortices in 2D 
hydrodynamics. 

It is well-known in astrophysics that a stellar system undergoes two successive types of relaxation (28j . In the 
"collisionless regime" , one can ignore correlations between stars due to finite N effects and make a mean field approx- 
imation. In that case, the evolution of the system is governed by the Vlasov equation. The Vlasov-Poisson system 
can undergo a violent collisionless relaxation towards a quasi stationary state (QSS) on a few dynamical times. This 
corresponds to the mixing of the distribution function by the mean field. This mixing process is responsible for 
the apparent regularity of galaxies. This is precisely what the Lynden-Bell [2^ statistical theory tries to describe. 
UnfortunatcUy, in astrophysics, the Lynden-Bell theory does not give a good prediction of the structure of the whole 
galaxy because of the problem of incomplete relaxation in the outer part of the system where mixing is inefficient. On 
a longer timescalc, in the so-called "collisional regime" , one must take into account correlations between stars due to 
finite N effects (graininess) (s^j- These correlations are expected to drive the system towards the ordinary statistical 
equilibrium state described by the Boltzmann distribution for t — > +00. This mixing process is due to discrete effects 
and takes place on a very long timescale. Indeed, the "collisional" relaxation time scales like {N/\nN)tjy and it di- 
verges for N +00 [2^. For self-gravitating systems, the Boltzmann distribution is not reached in practice because 
of evaporation [s^l and the gravothermal catastrophe [Sll, [s^] . Even if the statistical mechanics of self-gravitating 
systems is complicated because of the peculiar nature of the l/r interaction, the concepts of violent collisionless 
relaxation towards a QSS and of slow collisional relaxation towards the Boltzmann distribution are fundamental and 
can be transposed to other systems with long-range interactions. For example, they have been clearly illustrated for 
the HMF model [H. 

The same distinction between collisionless and collisional relaxation applies to the system of point vortices. In the 
"collisionless regime" , one can ignore correlations between point vortices due to finite N effects and make a mean 
field approximation. In that case, the evolution of the system is governed by the 2D Euler equation. The 2D Euler- 
Poisson system can undergo a process of violent collisionless relaxation towards a quasi stationary state (QSS) on a 
few dynamical times. This corresponds to the mixing of the vorticity by the mean field. This mixing process leads to 
the formation of large scale vortices similar to those seen in planetary atmospheres. This is precisely what the MRS 
statistical theory tries to describe. On a longer timescale, called the "collisional regime" , one must take into account 
correlations between point vortices due to finite N effects (graininess) [s^. These correlations are expected to drive 
the system towards the ordinary statistical equilibrium state, described by the Boltzmann distribution, for t +00. 
This mixing process is due to discrete effects and takes place on a very long timescale. As we shall see, the scaling 
with N of the "collisional" relaxation time of point vortices is not firmly established. 

It is of paramount importance not to confuse the Lynden-Bell (or MRS) statistical theory and the Boltzmann 
statistical theory which apply to very different timescales in the evolution of a Hamiltonian system with long-range 
interactions. This is related to the non-commutation of the limits N +00 and t — -l-oo. Mathematically speaking, 
the QSS is reached when the N +00 limit is taken before the t — ^ +00 limit while the Boltzmann distribution 
is reached when the t +00 limit is taken before the TV +00 limit. The Lynden-Bell (or MRS) theory makes 
the statistical mechanics of a continuous field evolving according to the Vlasov (or 2D Euler) equation while the 
Boltzmann theory makes the statistical mechanics of a discrete system of particles (stars or point vortices) evolving 
according to Hamiltonian equations. The first describes the QSS that is formed after a few dynamical times and 
the second describes the statistical equilibrium state that is reached for t — > +cxd. The distribution predicted by 
Lyndcn-BcU (or MRS) is different from the Boltzmann distribution due to additional conservation laws in the Vlasov 
(or Euler) equation. In fact, the Lynden-Bell (or MRS) statistical mechanics and the Boltzmann statistical mechanics 
describe two types of mixing occurring at different scales: in the process of violent "collisionless" relaxation, the 
smooth distribution function /(r,v,t) (or vorticity field uj{r,t)) that evolves according to the Vlasov (or Euler) 
equation is mixed by the mean field and the coarse-grained distribution function /(r,v,t) (or vorticity field aJ(r, t)) 
reaches a statistical equilibrium state. In the process of "collisional" relaxation, the discrete distribution of particles 
/d(r, V, t) = J2i fnS{r — ri(i))(5(v — v,;(<)) or a;d(r, t) = 'y6{r — ri{t)) (stars or point vortices) that evolves according 
to the Klimontovich equation (equivalent to the Hamilton equations) is mixed by discrete effects and the smooth 
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distribution of particles f{r,v,t) or uj{r,t) reaches a statistical equilibrium state. 

Onsager [l[ assumed ergodicity of the point vortex gas and determined the ordinary statistical equilibrium state 
by evaluating the density of states and the entropy of the Hamiltonian system. However, he did not mention the 
source of mixing leading to the statistical equilibrium state. As we have explained, there are two sources of mixing 
in the point vortex gas: the fluctuations of the mean field during the phase of violent relaxation and the fluctuations 
due to discrete effects during the phase of coUisional relaxation. The first one is very rapid and the second one very 
slow. It is likely that Onsager had in mind the mixing due to the fluctuations of the mean field. Indeed, he was 
basically interested in the 2D Euler equation, not by finite N effects ("collisions") in the point vortex gas. As we 
have explained, only continuous vorticity fields make sense in 2D hydrodynamics. The fact that he used the point 
vortex model was just a question of commodity. At that time, nobody knew how to make the statistical mechanics of 
a continuous vorticity fiow. However, the statistical equilibrium state of the 2D Euler equation should be described 
by the MRS distribution, not by the Boltzmann distribution. The Boltzmann distribution describes the statistical 
equilibrium state of the point vortex gas that is reached for t +00 due to the development of correlations (finite 
N effects) between point vortices. It is not clear whether Onsager had anticipated these two very different regimes 
related to the subtle non-commutation of the t +00 and N +00 limits. 

In the present paper, we shall be essentially interested in the slow "coUisional" relaxation of the point vortex gas 
towards the Boltzmann distribution due to finite N effects. Its fast "colhsionless" relaxation towards a QSS for 
N +00 will be only briefly discussed. Now, the relevance of the Boltzmann distribution for point vortices when 
t — > +00 is not clearly established because it relies on an assumption of ergodicity (or sufficient mixing) that may not 
be realized [s^ . If we want to prove that the point vortex gas truly relaxes towards the Boltzmann distribution for 
t +00, and if we want to determine the relaxation time (in particular its scaling with N), we must develop a kinetic 
theory. The MRS statistical theory is also based on an assumption of ergodicity. or efficient mixing, and the process 
of violent relaxation of the 2D Euler equation towards a QSS should also be described in terms of a kinetic theory 
(for the coarse-grained vorticity field) . This kinetic theory is more complicated to develop and will not be considered 
here. We only refer t o Issl-isst for some investigations in this direction. 

In previous papers |39l - l42j | . we have developed a kinetic theory of 2D point vortices and we have derived a kinetic 
equation for the evolution of the smooth vorticity field taking two-body correlations into account. This equation is 
valid at the order 1/iV and provides therefore the first order correction to the Euler equation obtained for N +00. 
This kinetic equation was derived by using various methods such as the projection operator formalism, the quasilinear 
theory and the BBGKY hierarchy. In these works, we focused on (distant) two-body collisions and neglected collective 
effects. This leads to a kinetic equation similar to the Landau [43 equation in plasma physics [qO]. A kinetic theory 
of the point vortex gas had been previously developed by Dubin & O'Neil [3] in the context of non-neutral plasmas 
under a strong magnetic field. They used the Klimontovich formalism and took collective effects into account. This 
leads to a kinetic equation similar to the Lenard-Balescu [45ll46| equation in plasma physics. Their work was pursued 
in (irl - fsot in different directions. In a recent paper [5l|, we used the Klimontovich formalism to derive a Fokker-Planck 
equation describing the evolution of a test vortex in a bath of field vortices, taking collective effects into account. 

The purpose of the present paper is two-fold. One, physical, motivation is to discuss the implications of the 
kinetic theory of two-dimensional point vortices. An important conclusion of our analysis is that the relevance of 
the Boltzmann distribution is not established by the present-day kinetic theory. More precisely, we show that, if the 
point vortex gas ever relaxes towards the statistical equilibrium state, this takes place on a very long time, larger 
than Nto (for axisymmetric distributions). The coUisional relaxation of point vortices is therefore a very slow process 
and requires to take into account three-body, four-body... correlations, of order that are ignored so 

far in the kinetic theory. It is also possible that mixing by "collisions" is not sufficient to drive the system towards 
the Boltzmann distribution. We emphasize, however, that the point vortex gas can rapidly reach a QSS as a result of 
a "colhsionless" violent relaxation. However, this QSS is described by the Miller-Robcrt-Sommeria, or Lyndcn-BcU, 
distribution, not by the Boltzmann distribution. Another, more technical, motivation of our paper is to derive the 
Lenard-Balescu-type kinetic equation for point vortices from the BBGKY hierarchy. This derivation (which constitutes 
the main result of the paper) is new and completes our previous derivation [4l| where collective effects were neglected. 
Our approach closely follows the method of Ichimaru |52| for deriving the Lenard-Balescu equation in plasma physics. 
We also consider the relaxation of a test vortex in a bath of field vortices described by a Fokker-Planck equation. 
Finally, we mention open problems and future directions of research. 

We note that other kinetic theories of point vortices exist but they apply to situations different from the one we will 
consider here. For example, Nazarenko & Zakharov [s^ obtained a kinetic equation for point vortices with different 
intensities moving on the background of a shear flow and experiencing "hard" collisions. Marmanis [s^ l and Newton 
& Mezic [sst derived kinetic equations for a vortex gas viewed as a coupling, via the Liouvillc equation, between 
monopoles, dipoles and tripoles. Finally, Sire & Chavanis [H^ developed a kinetic theory of three-body collisions 
(dipoles hitting monopoles) with application to the context of 2D decaying turbulence. General results about point 
vortices can be found in the book of Newton [53 . 
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II. TWO-DIMENSIONAL POINT VORTICES: EVOLUTION OF THE SYSTEM AS A WHOLE 

A. The Boltzmann distribution 

We consider an isolated system of point vortices with identical circulation 7 on an infinite plane. Their dynamics 
is fully described by the KirchhofF-Hamilton equations [s^, : 

dxi dH dyi dH -f'^ ^ 

^^=%' ^^ = -a^' H = --J2ln\r,^r,\, (1) 

i<j 

where the coordinates {x, y) of the point vortices are canonically conjugate. We shall denote the potential of interaction 
by u{r,r') ~ — ^ln|r — r'|. This Hamiltonian system conserves the energy E ~ H, the circulation F = Nj, the 
angular momentum L = ^ ■ jrf and the impulse P = J2i l^i ■ We take the origin of the system of coordinates at 
the center of vorticity so that P = (we shall ignore this constraint in the following). The angular momentum 
fixes the vortex size R. For t — >■ +00, we expect this system to reach a statistical equilibrium state described by the 
microcanonical distribution [l|, 0, |4l| : 

Pw(ri, rjv) = -7^J{E - H{r,, rr,))d{L - ^ jr^), (2) 

where g(E^ i) = / &{E — H)S{L — J2i 7^1) Yii is the density of states. The microcanonical entropy of the system 
is defined by S{E,L) = \ng{E, L). The microcanonical temperature and the angular velocity are then given by 
P = 1/T = {dS/dE)L and = {'^1 P){dS / dL)E- As first realized by Onsager the temperature of the point 
vortex gas can be negative. At negative temperatures (i{E) < 0, corresponding to high energy states, point vortices of 
the same sign group themselves in "supervortices" similar to large-scale vortices in planetary atmospheres. We define 
the thermodynamic limit as A'^ — ?> +00 in such a way that the normalized energy e = i?/F^, the normalized temperature 
77 = /3F7, the normalized angular momentum A = L/TR^ and the normalized angular velocity v = fl^R'^ /T are of 
order unity (these scalings result from simple dimensional analysis). We can renormalize the parameters so that the 
circulation of the vortices scales like 7 ~ and the vortex size like R ^ 1. Then, the inverse temperature scales 
like f3 ^ N, the energy like E 1, the angular momentum like L ^ I and the angular velocity like fi^ 1. On the 
other hand, the total circulation F ^ 1 and the dynamical time tjj ^ l/u R^/T ~ 1 are of order unity. 

In the thermodynamic limit A^ +00 defined previously, it can be rigorously proven ^^s^t the A^-body 

distribution function factorizes in a product of A^ one-body distribution functions 

TV 

PAr(ri,...,rw) = n^i(rO- (3) 

1=1 

Therefore, the mean field approximation becomes exact in the thermodynamic limit N — +00. Furthermore, the 
one-body distribution function Pi{r), or equivalently the smooth density of point vortices n(r) = NPi{r) or the 
smooth vorticity field w(r) = N'-fPi{r), is the solution of a maximum entropy principle 0, H^: 

S{E,r,L) ^max {S'sM | E[u}] = E , T[ij] = T , L[lj] = L} , (4) 

where 



Sb[lo] = - J -In-dr, E^-Ju^pdr, T = jujdY, L = J ujr^ dr, (5) 

are the Boltzmann entropy, the mean field energy, the circulation and the angular momentum. Here, is the mean 
field stream function produced by the smooth distribution of vortices uj according to the Poisson equation 

- AV^ = CO. (6) 

The mean velocity of a point vortex is then (V) = — z x Vip where z is a unit vector normal to the flow. Fundamentally, 
the Boltzmann entropy is defined by 5*5 = In where W is the number of microstates (complexions), specified by 
the precise position {ri, ...,rjv} of each point vortex, corresponding to a given macrostate, specified by the smooth 
vorticity field {a;(r)} giving the average number of point vortices in macrocells of size < A <^C 1. Using the Stirling 
formula for N ^ 1, wc obtain the expression ^ of the Boltzmann entropy. Introducing Lagrange multipliers and 



5 



writing the variational principle in the form 5Sb ~ I3SE — aST — I3^SL = 0, we obtain the mean field Boltzmann 
distribution 

^^^^e-/37(^+%.^)^ (7) 

where A is a positive constant. Substituting this relation in the Poisson equation ([6|), we obtain the Boltzmann-Poisson 
equation 

- AV' = ^76-^^^'^+^'''), (8) 

like in the theory of electrolytes (for /? > 0) or in the statistics of stars (for /3 < 0). The statistical equilibrium state 
is then obtained by solving this equation and relating the Lagrange multipliers /3, a (or A) and fi^ to the constraints 
E, r and L Then, we have to make sure that the resulting distribution is a maximum of Sb at fixed circulation, 
energy and angular momentum (most probable state), not a minimum or a saddle point. 



B. BBGKY-like hierarchy and 1/N expansion 

In order to establish whether the point vortex gas will reach the Boltzmann distribution ([7]) predicted by statistical 
mechanics and determine the timescale of the relaxation, in particular its scaling with N, we need to develop a kinetic 
theory of point vortices. Basically, the evolution of the A^-body distribution Ppf{ri, ...,rjv,i) of the point vortex gas 
is governed by the Liouville equation 

where 

^--^S'-'-i-ER^-i:^*-'-* (10) 

is the velocity of point vortex i due to its interaction with the other vortices. Here, tp{ri) denotes the exact stream 
function in produced by the discrete distribution of point vortices and V(j — > i) denotes the exact velocity induced 
by point vortex j on point vortex i. The Liouville equation ([3]), which is equivalent to the Hamiltonian system ([1]), 
contains too much information to be of practical use. In general, we are only interested in the one-body or two-body 
distributions [9l|. From the Liouville equation (|9]) we can construct a c omp lete BBGKY-like hierarchy for the reduced 
distributions Pj(ri, ...,rj,<) = / PAr(ri, r^, t) rfrj+i ...drAr. It reads [4]| : 

'-i+t t V{k^^)■^ + iN-J)±^ViJ + l^^)■^dr,^,=0. (11) 

i=l k=l,k=ii ' i=l * 

This hierarchy of equations is not closed since the equation for the one-body distribution Pi (ji , t) involves the two-body 
distribution P2(ri,r2,t), the equation for the two-body distribution P2(ri,r2,t) involves the three-body distribution 
P3(ri, r2, ra, t), and so on... The idea is to close the hierarchy by using a systematic expansion of the solutions in 
powers of 1/A^ in the thermodynamic limit N — ^ -l-oo. 
The first two equations of the hierarchy are 

BP f r)P 

i^(ri,t) + (7V-l) J V(2^1).^(ri,r2,t)dr2 = 0. (12) 

1 flP r)P r r)P 

-^(ri,r2,0 + V(2 ^ 1) . ^(ri,r2,t) + {N - 2) V(3 ^ 1) • ^(n, r2, rg, i) drg + (1 o 2) = 0. (13) 

2 Ot OTi J OTi 

We decompose the two- and three-body distributions in the suggestive form 

P2(ri,r2,t) = Pi(ri,t)Pi(r2,i) + P:;(ri,r2,i), (14) 



Ps (r 1 , r2 , rs , t ) - Pi (r 1 , t ) Pi (r2 , t )Pi (ra , t ) + P^ (r 1 , r2 , i )Pi (ra , t) 
+P;^ (r 1 , rg , t )Pi (r2 , t ) + P2' (r2 , rs , t )Pi (r 1 , t ) + P^ (r 1 , r2 , rg , t) . 



(15) 
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This is similar to the Mayer [s^ decomposition in plasma physics. Substituting Eqs. (|T4|) and (|T5|) in Eqs. p2|) and 
(|13p and simplifying some terms, we obtain 



9i 



(ri,t) + (7V-l) 



V(2^ l)Pi(r2,t)dr2 



9ri 



(ri,t) = -(iV-l) 



d 
dri 



V(2^1)P^(ri,r2,t)dr2, (16) 



1 9P' If 1 dP' 

-^(ri,r2,t) + (7V-2) y V(3^1)Pi(r3,t)dr3 •^(ri,r2,t) 

V(2^1)- ^V(3^1)Pi(r3,t)dr 

V(3^1)P^(r2,r3,t)dr3 



9Pi 

+ (A^-2 



^(ri,t)Pi(r2,t) + V(2 1) • ^(ri,r2,i) 

OTi OTi 



J V(3^1)P^(ri,r3,t)Pi(r2,t)dr3 + (iV-2)A. ^ V(3 ^ l)P^(ri, r2, r3, t) dr3 + (1 o 2) = 0. (17) 
and ([T71) arc exact for all N but they are not closed. We shall close these equations by expanding 



d_ 

dri 

Equations and ([T7|) are exact for all N but they arc not closed. We shall close these equations by expanding 
the solutions in powers of 1/A^ for ^ 1. In this limit, the correlation functions Pj(ri, r^, t) scale like 1/N^~^. 
In particular, P^ ~ 1/iV and Pg ~ 1/N'^. On the other hand. Pi ~ 1 and |V(i ^ j)| 7 - 1/iV. We are aimmg 
at obtaining a kinetic equation that is valid at the order 1/A^. Let us consider the terms in Eq. (|17p one by one. 
The first and second terms are of order 1/iV. They represent the transport of the two-body correlation function by 
the mean flow. The third term represents the effect of "soft" binary collisions between vortices; it is of order 1/iV. 
If we consider only these first three terms (as done in our previous paper [4l|), we obtain a kinetic equation that is 
the counterpart of the Landau equation in plasma physics. The fourth term represents the effect of "hard" binary 
collisions between vortices. This is the term which, together with the previous ones, gives rise to the Boltzmann 
equation in the theory of gases. It is of order l/N"^ but it may become large at small scales so its effect is not entirely 
negligible. For example, in plasma physics, hard collisions must be taken into account in order to regularize the 
logarithmic divergence that appears at small scales in the Landau and Lenard-Balescu equations. In the case of point 
vortices, there is no divergence at small scales in the kinetic equation that we shall obtain. Therefore, in this paper, we 
shall ignore the contribution of this term (but we note that it would be interesting to study it specifically). The fifth 
term is of order 1/iV and it corresponds to collective effects. In plasma physics, this term leads to the Lenard-Balescu 
equation. It takes into account dynamical Debye screening and regularizes the divergence at large scales that appears 
in the Landau equation. The main contribution of this work will be to take this term into account in the kinetic 
theory of point vortices in order to obtain a Lenard-Balescu-type kinetic equation from the BBGKY hierarchy. The 
last two terms are of the order 1/iV^ and they will be neglected. In particular, at the order 1/A^, we can neglect the 
three-body correlation function. In this way, the hierarchy of equations is closed. 

If we introduce the smooth vorticity field uj{ri,t) = N^Pi{ri,t) and the two-body correlation function ^(ri, r2,t) = 
iV^P2(ri,r2,t), we get at the order 1/A^: 



duj 
'dt 



(ri,i) 



N - 1 
N 



(V)(ri,t)-|^(ri,t) = -7 



-■ J V(2^1)g(ri,r2,t)dr2, 



(18) 



^||(ri,r2,t) + (V)(ri,t) ■ ^(ri,r2,t) 



i J V(3^1).g(r2,r3,t)dr3 



dri 



-^V(2 ^ 1) • ^cu{ri,t)Loir2,t) + (1 o 2) = 0. 

7"^ OTi 



We have introduced the mean velocity in ri created by all the vortices 

(V)(ri,t) = V(2 ^ I)c^(r2,t)dr2 = -z x ViP{ri,t), 
and the fluctuating velocity created by point vortex 2 on point vortex 1: 

V(2^1)=V(2^1)-l(V)(ri,0. 

We also recall that the exact velocity created by point vortex 2 on point vortex 1 can be written 

dui2 



V(2 -> 1) = -7z X 



dri 



(19) 



(20) 



(21) 



(22) 
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where U12 = w(|ri — r2|) is the binary potential of interaction between point vortices. Equations (|18|) and ([T9)) are 
exact at the order 0{1/N). They form the right basis to develop a kinetic theory of point vortices at this order of 
approximation . 



C. The limit A*' +00: the 2D Euler equation (collisionless regime) 

In the limit N — !■ +00 for a fixed time t, the correlations between point vortices can be neglected and the TV-body 
distribution factorizes in a product of N one-body distributions: 

N 

PN{ri,...,rN,t)^Y[Pi{r^,t). (23) 

Therefore, for long-range interactions, the mean field approximation is exact at the thermodynamic limit N — > +00. 
Substituting the factorization in the Liouville equation ©, and integrating on r2, ra, rjv we find that the 
smooth vorticity field uj{r, t) of the point vortex gas is solution of the 2D Eulcr equation 

^(ri,t) + (V)(ri,t).^(ri,t)==0. (24) 

This equation also results from Eq. ^TE\\ if we neglect the correlation function ^(ri, r2, t) in the right-hand side. The 
2D Euler equation describes the collisionless evolution of the point vortex gas for times smaller than Ntn- In practice, 
N is large so that the domain of validity of the 2D Eulcr equation is huge. The 2D Euler equation is the counterpart of 
the Vlasov equation in plasma physics and stellar dynamics. It can undergo a process of mixing and violent relaxation 
towards a quasistationary state (QSS) on a very short timescale, of the order of a few dynamical times to- This QSS 
has the form of a large-scale vortex. Miller and Robert & Sommeria [l^ have developed a statistical mechanics 
of the 2D Euler equation to predict these QSSs. The MRS theory is the counterpart of the Lynden-Bell theory for 
collisionless stellar systems [25|. The analogy between two-dimensional vortices and stellar systems is developed in 

M. 



D. The order 0{1/N): an exact kinetic equation (collisional regime) 



If we want to describe the collisional evolution of the point vortex gas, we need to consider finite N effects. Equations 
(fTS]) and HH) are valid at the order 1/7V so they describe the evolution of the system on a timescale of order Nto- 
The equation for the evolution of the smooth vorticity field is of the form 



duj, , - 



i(V)(ri,t).^(ri,t) = C[c.(ri,t)], 



(25) 



where C(u)) is a "collision" term analogous to the one arising in the Boltzmann equation in the theory of gases. In 
the present context, there arc no real collisions between point vortices. The term on the right-hand side of Eq. (|18p is 
due to the development of correlations between vortices. It is induced by the two-body correlation function .g(ri, r2, t) 
which is determined in terms of the vorticity by Eq. (fT9|) . Our aim is to derive a kinetic equation that is valid at the 
order 1/iV and that gives the first correction to the Euler equation. 
The formal solution to Eq. p9)) is 



(26) 



g(ri,r2,t) = -;^ J dr[ J dr'^ dt'Uiri,r[,t ~ t')Uir2,r'^,t - t') 



V(2'^l')-^+V(l' 



where the propagator U{ri, r[,t — t') satisfies the equation 
^ (r 1 , r; , i - i') + (V) (ri ,t) ■ ^U{r,,r[,t - t') ^ 



2') 



V(2^ l)C/(r2,r;,t-t')dr2 



(27) 



with the initial condition J7(ri, r'j^, 0) = S(vx — r'l). Equation (|27p can be viewed as a linearized version of the 2D 
Eulcr equation (sec Appendix IXj) . Indeed, if we make the replacement a; — !■ a; -f in Eq. (|A1[) and linearize it with 
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respect to 5uj [see Eq. (jA2p ]. we obtain Eq. 
obeys the linearized Euler equation. 

Substituting Eq. (|26|) in Eq. (jlSp . we obtain a kinetic equation 

d 



27|) in which U plays the role of 5uj. Consequently, the propagator U 



duj , 



^(V)(r„0.-(r,t)^-- 



dr- 



xV(2 



1) 



V(2' ^ 1') 



_9_ 



dr[ 
V(l' - 



dr'2 / drt/(ri,r'i,T)C/(r2,r'2,r) 



2') 



/o 



- T)w(r^,t - r), 



(28) 



that is exact at the order If we neglect collective effects, we recover the generalized Landau equation (|C6[) 

that was derived in our previous articles j4(]| - |42| . Equation (pS)) is a complicated non-Markovian integrodifferential 
equation. It is furthermore coupled to Eq. (|27p which determines the evolution of the propagator. In order to 
resolve this coupling, it is necessary to consider the timescales involved in the dynamics. We shall argue that, for 
a given vorticity profile w, the two-body correlation function relaxes to its asymptotic form on a timescale short 
compared with that on which to changes appreciably. This is the equivalent of the Bogoliubov ansatz in plasma 
physics. It is expected to be a very good approximation for iV ^ 1 since the two-body correlation function relaxes 
on a few dynamical times to while the vorticity field changes on a coUisional relaxation time of the order Nto or 
larger. Therefore, it is possible to neglect the time variation of uj{r,t — t) in the calculation of the collision term and 
extend the time integration to +00. This amounts to replacing the two-body correlation function in Eq. (jlSp by its 
asymptotic value (/(ri, r2, -I-cxd) for a given vorticity profile u. After the correlation function has been obtained as a 
functional of w, the time dependence of w can be reinserted. With this Bogoliubov ansatz (or adiabatic hypothesis), 
the kinetic equation (|28p can be rewritten 



'dt 



(ri,0 + 



N ~ 1 



TV 



-(V)(ri,0-|^(ri,0 



1 d 
1 dvi 



dv2 I dv'i I dv'2 



-t-00 



dTC/(ri,r'i,T)[/(r2,r^,T) 



xV(2 ^ 1) 



V(2'-l')-^+V(l'^2').A 



Similarly, the equation for the propagator takes the form 



97 



(ri 



(V)(ri,0-^C/(ri,r;,r) + 



V(2^ l)C/(r2,r'i,T)dr2 



dijj 



Lo{r[,t)uj{r'2,t). 



(ri,0 = 0, 



(29) 



(30) 



with the initial condition U{ri,r[,0) ~ d{ri — r[). The two equations (p9)) and (pO| are now completely decoupled. 
For a given vorticity profile a;(r, t) at time t, one can solve Eq. ([50)1 to obtain U{ri, r[, r) and determine the collision 
term in the right-hand side of Eq. (|29|) . Then, the vorticity profile Lu{r,t) evolves with time on a slow timescale 
according to Eq. (|29p . Interestingly, the structure of this kinetic equation bears a clear physical meaning in terms of 
generalized Kubo relations [4l[. This equation is valid at the order 1/TV and, for TV — > +00, it reduces to the (smooth) 
2D Euler equation ([24|) which describes the coUisionless evolution of the point vortex gas. 
The kinetic equation (j29p is, of course, equivalent to the pair of equations 



duj 
'dt 



(ri,t) 



TV - 1 
TV 



(V)(r„i).-(r„0 = -7^ 



V(2 ^ l)g(ri,r2,+oo)dr2. 



(31) 



ff(ri,r2, 



00) = - 
V(2' 



r 



dr[ / dr'^ 



i')-7|r + V(i'^2')-^ 



dT[/(ri,r'i,r)[/(r2,r'2,r) 
u;ir[,tMr'^,t). 



(32) 



which correspond to the first two equations of the BBGKY-like hierarchy at the order 1/TV within the Bogoliubov 
ansatz. These equations, supplemented by Eq. pop for the propagator, provide the formal solution of the problem in 
the general case. In order to obtain more explicit expressions, we have to consider particular types of flow. 



III. EXPLICIT KINETIC EQUATION FOR AXISYMMETRIC FLOWS 
A. Laplace-Fourier transforms 



We consider an axisymmetric distribution of point vortices that is a stable steady state of the 2D Euler equation. 
Therefore, the vorticity field evolves in time only because of the development of correlations between point vortices 
due to finite TV effects (graininess). In that case, an explicit form of the kinetic equation can be derived. 
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For an axisymmetric flow, introducing a system of polar coordinates, the vorticity field and the two-body correlation 
function can be written as uj{ri,t) = uj{ri,t) and g{ri,r2,t) = g{ri,r2,0i — 92,t), and the mean velocity as (V)(ri, t) = 
{V)g{ri,t)e0. On the other hand, according to Eq. ([22]) . the radial velocity (in the direction of ri) created by point 
vortex 2 on point vortex 1, is 

where ui2 = u(ri, r2,0i — O2) is symmetric in ri and r2 and even in (f) = 9i — 62 (see Appendix IbI) . In that case, Eqs. 
(PT|) and ^ take the form 




g{ri,r2,0i - 6*2, +00) 



i 1 r[dr[d9[J r'2dr'2de'2 dT^^{r[,r'2,9[ - e'2 



Id Id 



r'n dr'n 



7 ^{r[Mr'2) 
2 "'2/ J 

xU{n,r[,ei - 0[,T)U{r2,r'2, 62 - 0^, r). (35) 

For convenience, we have not written the time t in the vorticity field uj(r,t) appearing in the correlation function. 
As we have previously explained, the vorticity profile is assumed "frozen" on the short timescale that we consider to 
compute the asymptotic expression of the correlation function and the collision term (Bogoliubov ansatz). The time 
t will be restored at the end in the kinetic equation. 

We now expand the potential of interaction in Fourier series 

uir, r',e-e') = J2 e"(''-''^ii„(r, r'), (36) 

n 

and perform similar expansions for g{ri,r2,0i — 62) and U{ri,r'i^6i — di,t). In terms of these Fourier transforms, 
Eqs. ([M)) and ([55]) can be rewritten 

^(»"i,i) = 21717^ — — / r2dr2y^nUn{ri,r2)gn{ri,r2,+co), (37) 
dt ri dri Jq ^ 



gn{ri,r2,+oo) 



.(27r) 



7 Jo 



+ 00 



r[dr[ 



+ CXD 



+ CXD 



r2dr2 / dr nun{r'i^r'2) 



V,d7rV2d72)^^'^^^^'^\ 

xU„in,r[,T)U-nir2,r'2,T). (38) 



Introducing the Laplace transform of J7„ (ri , r[ , t) (see Appendix |A] for the definition of Laplace transforms) and 
integrating on time r, we get 



g«(ri,r2,-|-oo) 



-1 p+00 

- / r[dr[ / 
7 Jo Jo 



r'2dr'2 I da I da' - 



-nun{r^,r2) 



x[/„(ri,rl,a);7_„(r2,r^,a'), (39) 



where C is the Laplace contour in the complex cr-plane. The integration over a' can be performed by closing the 
contour by an infinite semicircle in the upper half-plane. Since U-n(r2,r'2,a') vanishes for \a'\ — > +cxd, the only 
contribution of the integral comes from the pole at cr' = —a. Using the residue theorem, we obtain 



/•+00 /•+00 /•+00 

5n(''i,''2,+oo) = / r[dr[ / rjdrj / danUn{r[ 

7 Jo Jo J-00 



' r' 
' '2 



xUn{ri,r[,a)U^n{r2,r'2, -<j). 



(40) 



Finally, substituting Eq. pO|) in Eq. (|37p . the kinetic equation takes the form 



duj 
'dt 



7-1 9ri Jo 



+00 



r-2 



(ir2 ^ nuniji, 



r2 



r[dr'i 



+00 



'2"' 2 



da nunir[,r'2, 



Uniri,r[,a)U-nir2,r'2, ~a). 



(41) 
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On the other hand, for an axisymmetric vorticity distribution uj(r), the Laplace-Fourier transform of the propagator 
is explicitly given by (see Appendix |^ : 

The first term on the right-hand side represents the advection by the mean flow, i.e. it corresponds to a pure rotation 
with angular velocity il{r). The second term takes into account collective effects. 

Before going further, some technical details must be given. Since [/„(ri,r^,a) is obtained as the Laplace transform 
of the propagator U„{ri, r[,T), this function is analytic in the upper half of the complex a plane. It is then continued 
analytically into the lower half-plane where it generally has singularities. The contour of a integration in Eq. (|4ip 
sees all singularities of [/^(ri, r'j^, cr) from above. On the other hand, the function C7_„(ri, r^, — cr) which is the complex 
conjugate to Un(ri,r'i, a) is an analytic function in the lower half-plane and is continued analytically into the upper 
half-plane where it generally has singularities. The integration contour in Eq. (|^T|) sees all the singularities of 
f7_„(ri,r2, — (t) from below. In order to take these boundary conditions into account, we shall write a + iO+ in place 
of a for the functions which are well defined and analytic in the upper half-plane and we use a — iO^ for the functions 
well defined and analytic in the lower half-plane. 

B. Without collective effects: Landau-type equation 

Before solving the general case, we first derive a Landau-type equation obtained by neglecting collective effects. In 
that case, the propagator (|42|) reduces to 

Unir,r\a) = ' (43) 

a — ml[r) Zirr 

and the kinetic equation (|4ip becomes 

dt ridriJo ^ (a - nS2(ri) zO+)((J - nr2(r2) - 20+) 

>^ l^TZ- - :^7Z-V(''i)^('^2). (44) 



T\ dri ?'2 dr 



The integration on a may be carried out by closing the contour with an infinite semicircle in the lower half-plane. 
Only the pole at cr = n^iri) contributes. Using the residue theorem, we obtain 



— (ri,t) = -27ri7— -— / r2 dr2 > n u„(ri,r2 
dt ri dri Jq ^ 



+°° — 2 1 ( \ d \ d 



/o ^ nQ,(ri) — nfl{r2) — iO+ \ri dri r2 dr2 

Then, with the aid of the Plemelj formula 



uj{ri)Lo{r2). (45) 



X ± iO+ V X 

the foregoing expression can be rewritten 



^ ='pf-\zfiT:S{x), (46) 



duj 1 d /■+°° /Id 1 d \ 

— {ri,t) = 2TT^^-— / r2dr2 y2n^Uniri,r2)^S[n{il{ri) - n{r2))] -7^ -r— w(ri)w(r2). 

dt ri dri Jo dri r2 dr2 } 



(47) 



Finally, using the identity (5(A.t) = ^(5(.t) and putting back the (slow) time dependence in the kinetic equation, we 
get 

dui 1 d r+°° f \ d \ d \ 

— (ri,t) = 2^2 — / r2dr2Y.\n\u^{ri,r2f5[n{rut)-^{r2.t)][—- —] Lo{ri,t)u{r2,t). (48) 

dt ri dri Jo „ "^i ^2 Of 2 / 

Using Eq. (jB4|) . the series can be explicitly calculated and we obtain the alternative form 

du) 1 d /■+°° /id 1 d \ 

— {ri,t) = 2TT^^-— / r2dr2xiri,r2)6[n{ri,t)-n{r2,t)]{-- \ oj{ri,t)Lo{r2,t), (49) 

dt ri dri J^ \ ri dri r2 dr2 J 
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with 



m— 1 ^ ^ 



i-(^ 



(50) 



where 



max(ri, and 



niin(ri,r2). This kinetic equation, which neglects cohective effects and takes into 



account only two-body encounters, is the counterpart of the Landau [43| equation in plasma physics. It was derived in 
[40l - |42j from different formalisms. In plasma physics, the Landau equation presents a logarithmic divergence at large 
scales. Therefore, collective effects which lead to Dcbye shielding, are crucial because they regularize the logarithmic 
divergence at large scales in the Landau equation. This is essentially what the works of Lenard [i^ and Balescu [i^ 
have demonstrated. Since the kinetic equation (|48)) does not present any divergence, the neglect of collective effects 
may not be crucial for point vortices. 



C. With collective effects: Lenard-Balescu-type equation 



If we take collective effects into account, the kinetic equation is obtained by substituting the expression (|^^ of the 
propagator in Eq. (|4ip and carrying out the integrations. The calculations are similar to those developed by Ichimaru 
[52| in his derivation of the Lenard-Balescu equation from the BBGKY hierarchy. 

Let us first perform the integration on r2- Using Eq. we obtain 



27r 



+ 00 



r2dr2 U„(ri,r2)f7_„(r2,r2, -cr) 



+ 00 



<5(r2-ry 



-(T + nQ{r2) r2 



- 2711 



G{-n,r2,r'2,-(T) 



{-a + nn{r2)){-a + nnir'^)) r2 dr2 



r2dr2 Un{ri,r2) 
I duj 

(r2) 



(51) 



which can be rewritten 



27r 



+ 00 



r2dr2 Un{ri,r2)U_n{r2,r'2, -cr) 



— (T + nfl{r!2) 



+ 00 



- / Mo/ ^ - / ^ G{-n,r2,r'^,-(T) duj ' 
Un[ri,r2) -211 / dr2Un[ri,r2) ; — —, — r — n—{r2) 



-a + nfl{r2) dr2 
Using the identity (jB8|) . this integral is finally reduced to the form 



2n 



r2dr2 u„(ri,r2)C/_„(r2,r2, -cr) = 



-cr + nfl(r'2) 



j--G{~-n,ri,r'2, -cr), 



(52) 



(53) 



where G{n, ri, r2, cr) is the "dressed" potential of interaction between point vortex 1 and point vortex 2. It is solution 
of the differential equation ()A13| . 

Substituting Eqs. ([^^ and ([5^ in the kinetic equation (PTj) . we obtain 



— (ri,i) = 27r7— — > 77, / r^dr^ / r2dr2 / danun[r-^,r2 

1 (5(ri - ri) G(n,ri,ri,cr) 



-n—^^{rx) 



cr-7if}(ri) + i0+ 27rri ((t - nf7(ri) + iO+)(cr - nf7(ri) + iO+) rxdrx 



12 



This equation can be decomposed as follows 



'dt 



1 d 
ri dri 



+ CXD 



+ CXD 



da nu„(ri,r2) 



Id 1 5 \ ■ 

— J. TTT '^(^1)^(^2) 

ri ori r'2 or'2 J 

G(-n,ri,r^, -a) 

(cr - nn{ri) + iO+)(cr - nl7(r^) - iO+) 



1 d 

ri ori 



+00 



da n 



{a - nVL{r-C} + iO+)(cr - nr2(r^) - iO+) 



w(ri)w(r^) 



r + 00 

x2tt I dr[un{r[,r'2)n—{r[) 



]__d_ 

ri dri 

G{n,ri,r[,a) 



-7— — Vn / r[dr[ / dan— 
n dri ^ 7o J -00 (o- 



G{n,ri,r[,a) 



nVL{ri) + iO+)(cr - nn{r[) + iO+) 

9r: 



— T^w(ri)w(r;) 



Using the identity (jBSp . we can write more succinctly 



1 d 
ri dri 



+ CX3 



+ 00 



dcr nun[ri,r'2) 



1 a Id 

ri dri 

G(-n,ri,r^, -cr) 



[a - nn{ri) + iO+)(cr - nf7(r^) - iO+) 



1 d 


ri dri 

1 d V- 

ri dri 



+00 



\ 2"' 2 



+00 



da 1 



G{-n,ri,r!2,-a) 



{a - nn{ri) + iO+){a - nr2(r^) - iO+) 



1 d 
ri dri 



oj{ri)uj{r'2) 



/ + 00 
da 1 
00 



G{n,ri,r[,a) 



X [-Un{ri,r'2) + G(n,ri,r2,CT)] 

(a - nO(ri) + iO+)(a - nn{r[) + iQ+) Vi^i'^^''^^'^^''^\ 
X [-u„(ri,ri) + G(-n,ri,r'i,-(T)] 



(55) 



(56) 



Parts of the first two terms cancel each other. On the other hand, the third term can be split into two parts. Hence, 
we get 



dijj 
'dt 



{rut) 



277 ri dri 



E^ 



da cj(ri) 



1 



a - nVL{ri) + iO+ 



„ 1^^°° , , , du , ,,G{-n,ri,r'-a) 



+7 



ri 9ri 

^_d_ 

ri dri 

1 a 



+ 00 



do n 



E^ 



+00 



r\dr\ 



— 00 
-\-QO 



da n 



1 u sr-^ 



+00 



r'ldr'i 



—00 
+00 



da 1 



G{~n, ri, rj, —a)G{n, ri, cr) 
(cr - nl7(ri) + iO+)(cr - nl7(r^) - iO+) 

G(n,ri,r^^,cr)Mn(r'^,ri) 

(cr - nf2(r-i) + iO+)(cr - nn{r[) + iO+) 

G{n,ri,r'i,a)G{~n,ri,r'i, -a) 
(cr - nn{ri) + iO+){a - nSl(ri) +iO+) 



cr - nf7(r;,) - tO^ 



i__d_ 

ri dri 



uj{ri)uj{r[) 



— ^uj{ri)uj{r'i) 
ri dri 



(57) 
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Using the identity (|B8[) . this can be rewritten 



duj 
'dt 



, . 7 1 9 ^-^ /■+°° , , , 1 
(ri,t) = —- — — / n da ulri) -—- 



' ri dri 

ri dri 

]__d_ 

ri dri 



+7— — X!*^ 



+°° G(-7i,ri,r^,-(T)G(n,ri,r^,cr) 
a(7 n 



+ 00 



r[dr[ 



da n 



+ 00 



r[dr[ 



— oo 

+ CJO 



da 1 



{a - nn{ri) + iO+){a - nf7(r^) - iO+) 

G{n,ri,r[,a)unir[,ri) 

(a - nO(r-i) + iO+){a ~ nn{r[) + iQ+) 

G{n,ri,r[,a)G{-n,ri,r[, -a) 
{a - nn{ri) + iO+){a - nn{r[) + iO+) 



(n) + iO+ 
X [-u„(ri,ri) + G{-n,ri,ri,-a)] 

ri ori 



']__d_ 

ri dri 

1 d 

ri dri 



uj{ri)uj{r'i) 
u}{ri)uj{r'i) 



(58) 



The integration in the first term may be carried out by closing the contour with an infinite semicircle in the lower 
half-plane; only the pole at cr = nfl(ri) — iO^ contributes. Similarly, by closing the contour in the upper half-plane, 
we find that the third term vanishes. The second and fourth terms may be combined with the aid of the Plemelj 
formula (^5)) . We are thus left with 



^(ri,i) = ij—-^y2nuj{ri) [-M„(ri,ri) + G(-n, ri, ri, -nfi(ri))] 

rii n Ori 



dt 



1 9 V- 2 /■ 

n dri ^ Jo 



^^^2 77^r~\ r,/ w I -04- ^w(ri w rz). 

n[il(r2j ~ it{ri)) + iO+ ri dri 



(59) 



Checking the symmetry of the terms with respect to the inversion of n, and again using the Plemelj formula (|46p. we 
get 



^(ri,t) = 7-^^^nw(ri)Im[G(n,ri,ri,nrj(ri))] 



ri dri ^ Jo 



1 d 



r2rfr2|G(n,ri,r2,nr!(r2))|^(5[n(f7(r2)-r!(ri))]-— w(ri)c^(r2) 

ri dri 



(60) 



where we have used G{—n,r,r', —a) = G*{n,r,r' ,a). To determine the first term in this equation, we multiply Eq. 
(|A13[) by 27rr G(— n, r, r', —a), integrate over r from to +oo and take the imaginary part of the resulting expression 
using the Plemelj formula (146]) . This yields 



ImG(n,r,r,cr) = -27r2 / dr' n\G{n,r,r' , a)\^S(a - 7in(r')) — (r'). 

Jo or' 



(61) 



Substituting this identity in Eq. (|60p and reinserting the (slow) time dependence in the the kinetic equation, we 
finally obtain 



duj 
'dt 



1 d r 

{ri,t) = 27r^^-— / r2dr2y2\n\\G{n,ri,r2,nniri,t))\^[niri,t) - nir2,t)] 
fi ori Jq 

Lu{ri,t)Lj{r2,t). (62) 



Id Id 

ri dri r2 dr2 



This can be rewritten 



did 1 d /■+°° 

— {ri,t) = 2TT^j-— / r2dr2xiri,r2,n{ri,t))S[n{ri,t)^n{r2,t)] 
dt ri dri Jq 



Id Id 
ri dri r2 dr2 



uj{ri,t)Lo{r2,t), (63) 



with 



X{'ri,r2) = ^ \n\\G{n,ri,r2,nn{ri,t))[' 



(64) 



14 



This kinetic equation, which properly takes collective effects into account, is the counterpart of the Lenard-Balescu 
equation in plasma physics. It was derived in [4^ [sij from the Klimontovich formalism. We have here provided an 
alternative derivation of this equation from the BBGKY hierarchy. Note that the Lenard-Balescu-type equation ([62|) 
differs from the Landau- type equation (j48p only by the replacement of the "bare" potential of interaction u„(ri,r2) 
by the "dressed" potential of interaction G(n, ri, r2, nf2(ri, t)) taking into account the contribution of the polarization 
cloud. This is similar to the case of plasma physics. 

Remark: the Landau-type kinetic equation (|48p and the Lcnard-Balcscu-type kinetic equation (|62p arc restricted to 
axisymmctric flows. By neglecting collective effects, we have obtained in [40l442| a generalized Landau equation (|C6[) 
that is valid for arbitrary flows. Sano (60j confirmed our results and attempted to derive a generalized Lenard-Balescu 
equation for arbitrary flows. Unfortunately, the final results are complicated and not fully explicit. Nevertheless, going 
beyond the assumption of axisymmetry as in [40l - |42j and [60| is certainly valuable in order to treat more realistic 
flows. The general kinetic equation ([^5]) could be a good starting point in this direction. 



D. Numerical resolution of the kinetic equation: No relaxation towards the Boltzmann distribution 

The kinetic equation ((63)) is valid at the order 1/A^ so it describes the "collisional" evolution of the point vortex 
gas on a timescale of order Nto- This kinetic equation conserves the circulation F, the energy E and the angular 
momentum L. It also monotonically increases the Boltzmann entropy in the sense that Sb > (_ff-theorem). These 
properties are proven in [4^ . l6l| . The change of the vorticity distribution in ri is due to a condition of resonance 
(encapsulated in the (5-function) between vortices located in ri and vortices located in r2 ^ ri which rotate with 
the same angular velocity D,{r2,t) = 0(ri,t) (the self-interaction at r2 = ri does not produce transport since the 
term in parenthesis vanishes identically) . Of course, this condition can be satisfied only when the profile of angular 
velocity is non-monotonic. The collisional evolution of the point vortices is thus truly due to long-range interactions 
since the current in ri is caused by "distant collisions" with vortices located in r2 7^ ri that can be far away. This is 
different from the case of plasma physics and stellar dynamics where the collisions are assumed to be local in space. 
The mean field Boltzmann distribution (O is a particular steady state of Eq. (|63p but it is not the only one: The 
kinetic equation ([63)) admits an infinite number of steady states. Indeed, all the vorticity profiles Lu{r) associated 
with a monotonic profile of angular velocity r2(r) are steady states of the kinetic equation (|63p since the (5- function is 
zero for these profiles. Therefore, the collisional evolution of the point vortex gas described by Eq. (|55)) stops when 
the profile of angular velocity becomes monotonic (so that there is no resonance) even if the system has not reached 
the Boltzmann distribution. In that case, the system settles on a QSS that is not the most mixed state predicted by 
statistical mechanics (sec Fig. [T]) [6l|. On the timescale Nto on which the kinetic theory is valid, the collisions tend 
to create a monotonic profile of angular velocity. Since the entropy increases monotonically, the vorticity profile tends 
to approach the Boltzmann distribution (the system becomes "more mixed" ) but does not attain it in general because 
of the absence of resonance. This is particularly obvious if we start from an initial condition with a monotonic profile 
of angular velocity that is non-Boltzmannian. In that case, the collision term vanishes, so that duj/dt = meaning 
that there is no evolution on the timescale Nt^. The Boltzmann distribution may be reached on longer timescales, 
larger than Nt^- To describe this regime, we need to determine terms of order N^"^ , or smaller, in the expansion of 
the solutions of the BBGKY hierarchy for N +00. This implies in particular the determination of the three-body, 
or higher, correlation function, which is a formidable task. At the moment, we can only conclude from the kinetic 
theory that, for an axisymmctric distribution of point vortices, the relaxation time satisfies 

tB. > NtD. (65) 

The collisional relaxation towards the Boltzmann distribution ([7]) is therefore a very slow process. In fact, up to 
now, there is no rigorous proof coming from the kinetic theory that the point vortex gas will ever relax towards 
the Boltzmann distribution predicted by statistical mechanics. Indeed, the point vortices may not mix well enough 
through the action of "collisions" . Therefore, the relaxation (or not) of the vorticity profile towards the Boltzmann 
distribution for t — -l-oo still remains an open problem. This is at variance with the Landau [43j and Lenard- 
Balescu [45I, llg equations of plasma physics which always converge towards the Boltzmann distribution (it is the 
unique steady state of these equations). Therefore, in plasma physics, it is sufficient to develop the kinetic theory 
at the order 1/iV. The kinetic theory of point vortices is more complicated since it requires to go to higher order 
in the expansion in power of 1/A^. This is similar to the case of spatially homogeneous one-dimensional systems 
with long-range interactions, such as one-dimensional plasmas [6^ IHl] and the HMF model [g^, [S^, for which the 
Lenard-Balescu collision term also vanishes at the order \/N . Therefore, tn > Ntu for these systems. For spatially 
homogeneous one-dimensional plasmas, the relaxation time is found numerically to scale like tji ~ N^tjj [66l.l67j which 
is the next order term in the expansion of the BBGKY hierarchy in powers of 1/N. For the spatially homogeneous 
HMF model, the scaling of the relaxation time with N is still controversial and different scalings such as ta ^ N'^to 
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FIG. 1: Evolution of the profile of angular velocity obtained by solving numerically the kinetic equation (|49|l [61l |. Time 
increases from bottom to top. The horizontal line shows the resonance between point vortices in r\ and which rotate with 
the same angular velocity. The evolution stops when the profile of angular velocity becomes monotonic so that there is no 
resonance anymore. The final profile does not correspond to the Boltzmann distribution ((T)). The kinetic equation (|49|l is valid 
on a timescale Nto- Therefore, on this timescale the vortex gas does not reach statistical equilibrium but remains "blocked" 
in a QSS with a monotonic profile of angular velocity. The Boltzmann distribution ((Tjl may be reached on a longer timescale 
(due to higher order correlations) but this timescale of relaxation is not known. It is not even clear whether the point vortex 
gas will ever relax towards the statistical equilibrium distribution ([7]) in the axisymmetric case. 



|6g | or tfi ^ e tjj [69| have been reported. An interesting problem would be to determine numerically the scaling 
with N of the relaxation time of an axisymmetric distribution of point vortices. If the collision term does not vanish 
at the next order of the l/7V-expansion, this would imply a timescale scaling like N'^tjj, but this scaling has to be 
ascertained (this project is currently under way). 

Note, finally, that the above results are only valid for axisymmetric distributions of point vortices. Non-axisymmetric 
distributions are described by a more complex kinetic equation (|29p which may present new resonances allowing the 
system to reach the Boltzmann distribution on a timescale of the order tn ^ Nt d (the natural first order of the kinetic 
theory). A linear Ntu scaling is indeed observed numerically for the relaxation of a non-axisymmetric distribution of 
point vortices [TOj . However, very little is known concerning the properties of Eq. (|29p and its convergence (or not) 
towards the Boltzmann distribution. It could approach the Boltzmann distribution (since entropy increases) without 
reaching it exactly. New resonances also appear for spatially inhomogeneous one dimensional systems with long-range 
interactions [4^ [tH . This may explain the linear Nt d scaling of the relaxation time observed numerically for spatially 
inhomogeneous one dimensional stellar systems (72| - [75j and for the spatially inhomogeneous HMF model [t^. On 
the other hand, for the HMF model, Yamaguchi et al. 73| find a relaxation time scaling like N^tjj with S ~ 1.7. 



In their simulation, the initial distribution function is spatially homogeneous but the coUisional evolution makes it 
Vlasov unstable so that it becomes spatially inhomogeneous. This corresponds to a dynamical phase transition from 
a non-magnetized to a magnetized state as theoretically studied in [t^. In that case, the relaxation time could be 
intermediate between N'^ta (permanently homogeneous) and Nto (permanently inhomogeneous). This argument 
(leading to 1 < (5 < 2) may provide a first step towards the explanation of the anomalous exponent S ~ 1.7 reported 
in (77} . The same phenomenon (loss of Euler stability due to "collisions" and dynamical phase transition from an 
axisymmetric distribution to a non-axisymmetric distribution) could happen for the point vortex system. 



IV. RELAXATION OF A TEST VORTEX IN A BATH: THE FOKKER-PLANCK EQUATION 



In the previous sections, we have studied the evolution of the system "as a whole" . In that approach, all the vortices 
are treated on the same footing. We now consider the relaxation of a "test" vortex (tagged particle) evolving in a 
steady distribution of "field" vortices. If the field vortices are at statistical equilibrium, described by the Boltzmann 
distribution (O, their density profile does not evolve at all. In that case, they form a thermal bath. However, we shall 
also consider the case of an out-of-equilibrium (i.e. non-Boltzmannian) bath corresponding to a vorticity distribution 
U!{r) that is a stable steady state of the 2D Euler equation with a monotonic profile of angular velocity. As we have 
seen previously, this distribution does not change on a timescale of order Nt^. Since the relaxation time of a test 
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vortex in a bath is of order {N/ In N)ti:) (see below) we can consider that the distribution of the field vortices is 
"frozen" on this timescale. They form therefore an out-of-equilibrium bath. 

Let us call P(r, t) the probability density of finding the test vortex at position r at time t. For simplicity, we 
shall consider axisymmetric distributions (some results valid in the general case can be found in [4l|). The evolution 
of P(r, t) can be obtained from the kinetic equation (1631) by considering that the distribution of the field vortices, 
described by the vorticity profile w(r2,i), is fixed [92|. In the BBGKY hierarchy, this amounts to specializing a 
particular point vortex in the system (the test vortex described by the variable 1) and assuming that the other 
vortices (the field vortices described by the running variable 2) are in a steady state. If we replace u){ri^t) by P(r,t) 
and uj{r2,t) by uj{r'), we get 



^{r,t) = 27t'i~ r'dr'x{r,r'Mr))6{n{r) ~ n{r')) - -^j P(r,i)w(r'). 



(66) 



This procedure transforms the integrodiffcrential equation describing the evolution of the system "as a whole" 
into a differential equation (|66|) describing the evolution of a test vortex in a bath of field vortices. Equation (|66|) can 
be written in the form of a Fokker-Planck equation 



dP 
'dt 



ld_ 

r dr 



involving a diffusion coefficient 



Dir) 



+ 00 



and a drift term due to the polarization 



r'dr'xir, r' ,n{r))6{n{r) - 17(r'))w(r'), 



+ 00 1 

dr'xir,r' ,n{r))S{nir) ~ nir')) — {r'). 

dr 



(67) 



(68) 



(69) 



Physically, the diffusion coefficient is due to the fiuctuations of the velocity field produced by a discrete number of 
point vortices; it can be directly derived from the Kubo formula [1^. On the other hand, the drift arises from the 
retroaction of the perturbation on the field vortices induced by the test vortex, just like in a polarization process; it 
can be directly derived from a linear response theory [39j . In the present case, the coefficients of diffusion and drift 
depend on the position r of the test vortex. Hence, it is more appropriate to write Eq. ()66p in a form which is fully 
consistent with the general Fokker-Planck equation 
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(Ar) 
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drift 



dD 
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(70) 



(71) 



Substituting Eqs. and in Eq. ([TTt . and using an integration by parts, we find that the total drift term is 
given by 



V;''"/*(r) = 2Tr^j / rr'dr'uj{r') 



+ 00 



ld_ 

r dr 



1 d 
r' dr' 



x{r,r',n{r))S{n{r)~n{r')) 



1 



(72) 



The two expressions (|57)) and ([701 of the Fokker-Planck equation have their own interest. The expression ([70)1 where 
the diffusion coefficient is placed after the second derivative d^ (DP) involves the total drift "Vdri ft and the expression 
(p7)) where the diffusion coefficient is placed between the derivatives dDdP isolates the part of the drift Vpoi due to 
the polarization. This expression is directly connected to the form of the Lenard-Balescu-type equation ((63l) . It has 
therefore a clear physical interpretation. The Fokker-Planck equation (j66p and the expressions (|68|) . (|69l) and ([72]) of 
the diffusion coefficient and drift term can also be obtained directly by calculating the first and second moments of 
the increment of radial position of the test vortex using the Klimontovich approach [5l| . 

If the profile of angular velocity of the field vortices ^{r) is monotonic, we can use the identity d{il{r) ~ il{r')) ~ 
S{r - r')/\n'{r)\ and we find that 



Dir) 



2 Xjr, r,n{r)) 
^ 



a;(r). 



(73) 
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and 



" ^ |E(r)| dr^^' 

where S(r) = rfl'(r) is the local shear. If we neglect collective effects, we can replace xi'^'i'^'i^i''')) by 



(74) 



1 1 1 

c{r,r) = V |7i|t/2 (r,r) = ^ - = ^liiA, 



(75) 



ri=l 



where InA = "^^=1 n ''^ Coulomb factor 93 1 that has to be regularized with appropriate cut-offs as discussed in 



[33, [so, IMI ■ It is then found that In A ^ ^ln^ in the thermodynamic limit N +oo. Wc note that the diffusion 
coefficient and the drift due to the polarization are inversely proportional to the shear. Furthermore, the diffusion 
coefficient is proportional the vorticity profile of the field vortices while the drift is proportional to its gradient. 
Comparing Eqs. ([75)1 and ([71)1 . we find that the drift velocity is related to the diffusion coefficient by the relation 



VP°\r) 



dlnuj 
D{r)——{r). 
ar 



(76) 



This can be viewed as a generalized form of Einstein relation for an out-of-equilibrium distribution of field vortices. 
Combining the previous results we find that the Fokker-Planck equation (|67p can be written 



dP 
'dt 



ld_ 

r dr 



-D{r) 



\ dr dr 



(77) 



with a diffusion coefficient given by Eq. ([75]) . 

If the field vortices are at statistical equilibrium (thermal bath), their vorticity profile is the Boltzmann distribution 



uj{r) ^ A-fe-^^'^'''^''\ 

where ^J^ir) = ip{r) + ^r^ is the relative stream function (see Sec. Ill Ap . We have 

^(r') = -^jcoir'Ar') = /?7^(r')(f7(r') - f^^)/, 
dr dr 



(78) 



(79) 



where we have used fl{r) = —{l/r)dip/dr. Substituting this relation in Eq. ((69)) . using the (5-function to replace ^{r') 
by r2(r), using il(r) — flL — —{\/r)d%l)^/dr and comparing the resulting expression with Eq. (|68p . we finally obtain 



Vr\r) = -D{r)P^'^{r). 



(80) 



The drift velocity Ypoi = — -D/37V?/'* is perpendicular to the relative mean field velocity (V,) = — z x V■^p^,. Further- 
more, the drift coefficient (or mobility) satisfies an Einstein relation ^(r) = D{r)l3^. We recall that the drift coefficient 
and the diffusion coefficient depend on the position r of the test vortex and that the temperature is negative in cases 
of physical interest. We also stress that the Einstein relation is valid for the drift 'Vpoi due to the polarization, not for 
the total drift W dri ft which has a more complicated expression. We do not have this subtlety for the usual Brownian 
motion where the diffusion coefficient is constant. For a thermal bath, using Eq. ([SO)) , the Fokker-Planck equation 
(p7)) can be written 



dP 
'dt 



ld_ 

r dr 



, JdP dtb^ 



(81) 



where D{r) is given by Eq. with Eq. ([75)1 . Of course, if the profile of angular velocity of the Boltzmann 

distribution is monotonic, we find that Eq. ([77)) with Eq. ([75[) returns Eq. ([ST]) with a diffusion coefficient given by 
Eq. (|73p with Eq. ([75[). Finally, we note that the systematic drift VpoZ = —DjS^Vip^, of a point vortex [S^] is the 
counterpart of the dynamical friction Fpoj = —D\\j3mw of a star [79| . Similarly, the Smoluchowski-type form of the 
Fokker-Planck equation (I5T]) describing the relaxation a point vortex in a "sea" of field vortices [39| is the counterpart 
of the Kramers- type form of the Fokker-Planck equation describing the relaxation of a star in a globular cluster 7^ . 
This is an aspect of the numerous analogies that exist between two-dimensional vortices and stellar systems (26| . 
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The Fokkcr-Planck equations ((77)) and (|8T|) have been studied for different types of bath distribution in [6l| . The 



distribution of the test vortex P(r, t) relaxes towards the distribution of the bath w(r)/r on a typical timescale 

tT' - ^tD, (82) 

where the logarithmic correction comes from the scaling with N of the Coulombian factor In A. However, the relaxation 
process towards the steady state is very peculiar and differs from the usual exponential relaxation of Brownian particles. 
In particular, the evolution of the front profile in the tail of the distribution is very slow, scaling like (Int)^/^, and the 
temporal correlation function (r(O)r(i)} decreases algebraically, like \nt/t (for a thermal bath), instead of exponentially. 
This is due to the rapid decay of the diffusion coefficient D{r) for large r. Similar results had been found earlier for 
the HMF model [sllH^. 

Finally, we stress that the evolution of the system "as a whole" is very different from the evolution of a test vortex 
in a bath. We have seen in Sec. IIIIDI that the relaxation time of the system as a whole is strictly larger than Ntjj 
(for axisymmetric distributions) while the relaxation time of a test vortex in a bath is of the order (N/ lTiN)tD- In 
particular, a steady state of the 2D Euler equation with a monotonic profile of angular velocity does not change 
on this timescale. This justifies our procedure of developing a bath approximation for out-of-equilibrium (i.e. non- 
Boltzmannian) distributions of the field vortices. 



V. CONCLUSION 

The statistical mechanics of two-dimensional point vortices was first considered by Onsager [l| in a seminal paper. 
The point vortex gas is a iV-body Hamiltonian system with long-range interactions that shares many analogies with 
electric charges in a plasma and stars in a stellar system. Like stellar systems, the point vortex gas displays two 
successive types of relaxation: A fast, or violent, relaxation due to mean field effects followed by a slow relaxation due 
to discrete effects. 

In a first regime, before the correlations ("collisions") due to graininess and finite N effects comes into play, the 
evolution of the smooth vorticity field is governed by the 2D Euler equation. In practice, the validity of the 2D 
Euler equation is huge since the coUisional relaxation time diverges with the number of point vortices. Starting 
from a generically unstable or unsteady initial condition, the 2D Euler- Poisson system experiences a process of violent 
relaxation leading to a non-Boltzmannian quasi-stationary state (QSS) which is a steady state of the 2D Euler equation 
on the coarse-grained scale. These QSSs correspond to large-scale vortices observed in 2D turbulence like monopoles, 
dipoles, tripoles... One can attempt to predict these QSSs in terms of statistical mechanics by using the Miller- 
Robert- Sommeria [l3, [13 theory which is the hydrodynamical version of the Lynden-Bell (2^ theory. This theory is, 
however, based on an assumption of ergodicity (or efficient mixing) which may not always be fulfilled. This leads to 
the complicated problem of incomplete relaxation which restricts the range of validity of the MRS (or Lynden-Bell) 
statistical theory. For example, "vortex crystals" [8l| could correspond to metastable states that are only partially 
mixed. 

On a longer timescale, "distant collisions" (i.e. correlations, graininess or finite N effects) between vortices must be 
taken into account. At the order 1/A^, this collisional evolution is described by the kinetic equation (|29p which reduces 
to the more explicit form (j63|) for an axisymmetric distribution of point vortices. However, this kinetic equation does 
not relax towards the Boltzmann distribution predicted by statistical mechanics for t — >■ -l-oo. Indeed, the evolution 
stops when the profile of angular velocity becomes monotonic. The relaxation towards statistical equilibrium may 
occur on a very long time, larger than Nto- To settle this issue, we have to develop the kinetic theory at the order 
1/iV^, 1/N^,... by taking into account three-body, four-body,... correlation functions. It is also possible that the point 
vortex gas never reaches the Boltzmann distribution: Mixing by "collisions" may not be efficient enough. For non- 
axisymmetric distributions of point vortices, the relaxation time may be reduced, due to the occurrence of additional 



resonances, and achieve the natural timescale Nto predicted by the first order kinetic theory [42|. This linear scaling 
has been observed in numerical simulations (TOj . It is, however, hard to prove that the general kinetic equation (|29p 
will reach the Boltzmann distribution. It may approach it without reaching it exactly. 

Finally, the stochastic motion of a test vortex in a "sea" (bath) of field vortices can be described in terms of a 
Fokker-Planck equation involving a diffusion term and a drift term. For a thermal bath, they are connected to each 
other by an Einstein relation. The diffusion coefficient (resp. drift term) is proportional to the vorticity distribution 
(resp. to the gradient of the vorticity distribution) of the field vortices and inversely proportional to the local shear. 
The distribution of the test vortex relaxes towards the distribution of the field vortices on a timescale of the order 
{N/ \nN)tD but the relaxation process is peculiar [Glf . 

For the point vortex gas, it is important to clearly distinguish the phase of violent collisionlcss (correlationless) 
relaxation towards a non-Boltzmannian QSS, taking place on a few dynamical times tu, from the slow collisional 
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(correlational) relaxation towards the Boltzmann distribution taking place on a very long timescale > Ntjj. In 
this sense, the MRS theory and the Boltzmann theory describe two completely different regimes. This distinction 
between coUisional (Boltzmann) and collisionless (Lynden-Bell) relaxation has been clearly understood for a long 
time in astrophysics [28l . [29| but it may be less well-known in the context of two-dimensional point vortex dynamics. 
This is because, implicitly, the Onsager theory has been used to describe the equilibrium state resulting from the 
phase of violent relaxation before the Lynden-Bell and the MRS theory were formulated. In fact, the Boltzmann 
distribution resulting from the Onsager theory has two completely different interpretations: (i) it can be viewed as 
an approximation of the MRS distribution arising from the collisionless mixing of the 2D Euler equation on a "short" 
timescale; (ii) it can be viewed as the ordinary Boltzmann equilibrium state arising from the coUisional evolution of 
the point vortex gas on a long timescale. In the first case, mixing is due to a purely mean field process while, in the 
second case, mixing is due to discrete interactions between vortices (collisions). These two regimes correspond to a 
different ordering of the N — +00 and t — )- +00 limits. 

We must also clearly distinguish the evolution of the system of point vortices "as a whole" from the relaxation 
of a test vortex in a bath. They correspond to completely different situations and they are described by different 
equations. The evolution of the system as a whole is described by the integrodifferential equation (|63| similar to the 
Lenard-Balescu equation in plasma physics, while the relaxation of a test vortex in a bath is described by a differential 
equation (j66[) similar to the Fokker-Planck equation in Brownian theory. 

Finally, we may compare the evolution of a large, but finite, number of point vortices with the evolution of a real 
continuous vorticity field. On a "short" timescale, we can neglect correlations and finite N effects in the evolution of 
point vortices and we can neglect viscosity in the evolution of the real continuous vorticity field. In that case, the two 
systems can be approximated by the 2D Euler equation and they should behave similarly. In particular, they may 
experience a process of violent relaxation towards a QSS. On a "long" timescale, the evolution of the point vortex 
gas and of the real continuous vorticity field will differ. Point vortices will feel correlations due to finite N effects and 
will (with the reserve that we have given in this paper) relax towards the Boltzmann distribution while a continuous 
vorticity field will feel inherent viscosity and will decay to zero. 



Appendix A: The linearized 2D Euler equation 



The 2D Euler equation can be written 



^ + u-Vw = 0, (Al) 
at 



where u = — z x Vip is the velocity field (z is a unit vector normal to the plane of the fiow). The stream function '0(r, t) 
is related to the vorticity u;(r, t) by the Poisson equation lS.1l} = —u. The potential of interaction it(r, r'), which is the 
Green function of the Laplacian operator A, is defined by Au(r,r') = —6{r — r'). In an infinite domain, u ~ w(|r — r'|) 
only depends on the absolute distance between two points and is given by u(|r — r'|) = — (l/27r)ln|r — r'|. This 
corresponds to a Newtonian (or Coulombian) interaction in two dimensions. Therefore, the stream function is related 
to the vorticity field by an expression of the form '(/'(r, t) = / M(|r — r'|)a;(r', t) dr. This can be written as a product 
of convolution ip ~ u * lo. 

Let us consider a small perturbation 6uj{r,t) around a steady state a;(r) of the 2D Euler equation. We write 
uj{r,t) = uj{r) + Suj{r,t), ip{r,t) = V'(r) -I- Stp{r,t) and u{r,t) = u(r) -I- Ju(r,t). Substituting these decompositions in 
Eq. (jAip and neglecting the quadratic terms, we obtain the linearized Euler equation 

— ^ + u • VSlj + Su • Vcj = 0. (A2) 
dt 

If we restrict ourselves to axisymmetric mean fiows, then u(r) = u{r)e0 with u(r) = ~^('') — ^i'r)r, where 
^{f) — -p- Jq u}{r')r' dr' is the angular velocity (see Sec. 6.1 of [6l|)- In that case, Eq. (jA2p becomes 

dSuj ^, .dSuj IdSipduj 

To solve the initial value problem, it is convenient to introduce Fourier-Laplace transforms. The Fourier-Laplace 
transform of the perturbation of the vorticity field duj is defined by 

/•2Tr in r+oo 

dQin,r,a)= — dte-'^''^-"'Hu{e,r,t). (A4) 

./n 27r In 
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This expression for the Laplace transform is vahd for lm{a) sufficiently large. For the remaining part of the complex 
cr plane, it is defined by an analytic continuation. The inverse transform is 



+00 „ , 



(A5) 



where the Laplace contour C in the complex a plane must pass above all poles of the integrand. Similar expressions 
hold for the perturbation of the stream function Sip. If we take the Fourier-Laplace transform of the linearized Euler 
equation (|A3[) . we obtain 



(A6) 



1 duj 

— Su)(n, r, 0) — iaSujin, r, a) + inQMSuin, r, a) + in——-6'ijj(n, r, a) = 0, 

r or 

where the first term is the spatial Fourier transform of the initial value 

/'27r in 

Suin, r, 0) = / — e~""'%u;{e, r, 0). 
./n 27r 



(A7) 



The foregoing equation can be rewritten 



, 1 duj 



ouj[n,r,aj = —j!---—d'ip(n,r,a) 



(7 — nil(r) 



i{(7 — nf2(r)) ' 



(A8) 



where the first term on the right hand side corresponds to "collective effects" and the second term is related to the 
initial condition. The perturbation of the stream function is related to the perturbation of the vorticity by the Poisson 
equation 



r dr dr 89^ 
Taking the Fourier-Laplace transform of this equation, we obtain 



(A9) 



1 d d v? 
r dr dr 



Substituting Eq. in Eq. ((XTOl) . we find that 



Sip{n, r, cr) = —5Lb{n, r, a). 



(AlO) 



\ d d n 



1 1 ^ 

' r dr 



r dr dr r^ a — nfl{r) 



Sip{n, r, a) 



5u!{n, r, 0) 
i{a — nQ{r)) 



(All) 



Therefore, the Fourier-Laplace transform of the perturbation of the stream function is related to the initial condition 

by 



S'ip{n,r,(T) = - I 2'Kr' dr' G{n, r, r',a) ^ , 



where the Green function G(ri, r, r', cr) is defined by 
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r dr dr r^ cr — nQ(r) 



G{n,r,r',a) = 
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If we neglect collective effects in the foregoing equation, we obtain 



1 d d n 
r dr dr r 



Gbare{n,r, r') 



5{r-r') 
2iTr 



(A14) 



Therefore, Gbare{n-,r,r') = u„(r, r') is the Fourier transform of the "bare" potential of interaction u that is solution 
of the Poisson equation Au = ~d, while G(n, r, r', cr) is the Laplace- Fourier transform of the potential of interaction 
"dressed" by its polarization cloud, i.e. taking collective effects into account. 
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Using the previous formulae, we can now relate the Laplace-Fourier transform <5(2;(n, r, a) of the perturbation of the 
vorticity to the initial condition Suj{n,r,0). Substituting Eq. (|A12p in Eq. (|A8p . we obtain 

If wc consider an initial condition of the form 5uj{r, 9, 0) ~ ^5(9 — 6')S{r — r')/r, implying 

6oj{n, r, 0) = p- '''''' S{r - r'), (A16) 
27rr 



we find that 



n 



1 d, 



5uj{n,r,a) = ^ -/c^ ^ v > > / „, (A17) 

' ' a — nfl(r) i(cr — nCl(r')) 27rr i(cr — nfl(r)) 

As shown in Sec. IIIIl this equation gives the expression of the propagator C/„(r, r', a) [see Eq. This operator is 

also called the resolvent operator as it connects Suj{n, r, a) to its initial value: 

r + co 

6u;{n,r,a)^ U„ir,r" ,a)6uj{n,r" ,0)2T:r" dr" . (A18) 



Appendix B: An integral equation 

The perturbation of the stream function is related to the perturbation of the vorticity by an expression of the form 

6i}{r,t)^ I u{\r -r'\)6uj{r',t)dr. (Bl) 



The potential of interaction can be written 

u(|r - r'l) = u (^y^r^ + r'^ - 2rr' cos{e - 6*')) = u{r, r' , <j)), (B2) 

where (j) = 9 — 9' . We note that u{r,r'(f>) ~ u{r\r,(j)) and u{r,r' , —(j)) = u{r,r',(l>). Due to its 0-periodicity, the 
potential of interaction can be decomposed in Fourier series as 

uir, r', 0) = V e™*w„(r, r'), u^{r, r') = ^ ^u{r, r' , </)) cos(n(A). (B3) 

For the Coulombian potential it(|r — r'|) = — (l/27r) In |r — r'|, satisfying u{r, r',(l>) = — (l/47r) ln(r^ + r'^ — 2rr' cos (/>), 
the integrals in Eq. (jB3[) can be performed analytically [i^ . We find that 

1 /r \ 1 
^in(r,r') = ^(^^j ' Wo(r,r') = -— lnr>. (B4) 

Therefore, the potential of interaction can be written 

|n| 



w(r,r',0) = --!-lnr> + j- 7^ ( — ) 



Jn<l> 



(B5) 



which is just the Fourier decomposition of the logarithm in two dimensions. Taking the Fourier-Laplace transform of 
Eq. (|Bip and using the fact that the integral is a product of convolution, we obtain 

p+oo 

S'ip{n,r,a) ^ 2tt / r'dr' Un{r,r')Suj{n,r' ,a). (B6) 



(B7) 



If we substitute Eq. (jA8[) in Eq. (|B6|) . we obtain the equation 



Stp(Ti,r,<T) — 2tt / r' dr' ■Un(r,r') — '' 6ip(n,r' ,a) — —2n / r'dr' Un(r,r')—— 

Jo cr - ni2(r') Jq H^^ 



nr2(r'))' 
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which is equivalent to Eq. (jA12[) . This imphcs that the Green function G{n,r,r' ,a) satisfies an integral equation of 
the form 



G{n, r, r',a)-2TT r" dr" u„ {r, r") 

Jo CF - nil[r") 



G{n,r",r',a) = u„(r,r'), 



which is equivalent to Eq. (jAlSp . If we neglect collective effects, Eq. (|B7[) reduces to 



Alt ^ 9 ^ ,jLo{n,r,0) 

oip[n, r, a) ~ — 27r / r dr Un(r, r )- 

Jo 



(B8) 



(B9) 



Comparing Eq. (|B9|) with Eq. (|A12p . we see that the bare Green function Gbarein.r.r') — w„(r, r') is the Fourier 
transform of the potential of interaction u. 



Appendix C: General kinetic equation without collective effects 

If we neglect collective effects, the first two equations and ([T^ of the BBGKY-like hierarchy reduce to 

|^(ri,t) + ^(V)(ri,t) . ^(ri,t) = -7^ ' / ^ l)ff(ri, r^, t) dr^, (CI) 



|(rr,r2,t) 



1 

'7 



V(2^1). A+v(1^2). A 

OTi OV2 



Lu{ri,t)uj{r2,t). 



Equation (jC2p is just a first order differential equation in time. It can be written symbolically as 

dg 



dt 



Cg = S[uj\, 



(C2) 



(C3) 



where £ is an advection operator and 5 is a source term: The correlation function is transported by the mean flow 
(l.h.s.) and modified by two-body collisions (r.h.s.). This equation can be solved by the method of characteristics. 
Introducing the Green function 



G(t,t') =exp|- jj" /:(r)rfrj 



constructed with the smooth velocity field (V), we obtain 



ff(ri,r2,i)== dTG(t,t-T) 

7 Jo 



■v(2-l)-^+V(1^2).A 



uj{ri,t - T)w(r2,t - t), 



(C4) 



(C5) 



where we have assumed that no correlation is present initially so that g(ri,r2,t = 0) = (if correlations are present 
initially, it can be shown that they are rapidly washed out). Substituting this result in Eq. (jCip . we obtain 



du! 
'dt 



^(V)(r„t).— (r„^).- 



^ / dr /" dr2V'^{2^l)G{t,t~T) 
1 "'0 J 



^"^(2^ 1)71^ + ^^(1 ^2) A 



Uj{ri,t- T)-{V2,t- t). 

1 



(C6) 



In writing this equation, we have adopted a Lagrangian point of view: the coordinates following the Greenian must 
be viewed as Ti{t — t)^ Yi{t) — JJ" ds (V){ri{t — s),t — s) ds. The generalized Landau equation (jC6p which is valid for 
possibly non-axisymmetric flows and which takes into account non-Markovian effects has been derived in our previous 
papers from the projection operator formalism [l^l, the BBGKY hierarchy [4l| and the quasilinear theory based on 
the Klimontovich equation |41| . The generalized Lenard-Balescu equation (|28|) can be viewed as an extension of this 
equation taking collective effects into account. Note that the Markovian approximation may not be justified in every 
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situation since it has been found numerically that point vortices can exhibit long jumps (Levy flights) and strong 
correlations js^ [ssj. 

If we implement the Bogoliubov ansatz and extend the integration on time r to infinity, we get 



Q r+oo . 



a;(ri,t)-(r2,t), 

7 



(C7) 



where, no w, r,; (t — r) = ri{t) — JJ" o?s ( V) (r^ (<: — s),t)ds. If we consider axisymmetric flows, this equation can be 
simplified j40l - l42l | and we obtain the Landau-type equation If we take collective effects into account, Eq. (|C7[) 

is replaced by Eq. ((29)) which reduces to the Lenard-Balescu-type kinetic equation (|62|) for axisymmetric flows. 
If we assume right from the beginning that the mean flow is axisymmetric, Eqs. (|CI|) and (|C2p can be rewritten as 



du! 
'dt 



{rut) 



ri dri 



+ 00 



r2 dr2 



d02 T^(ri,r2, 







(C8) 



If (ri , r2 , 01 - 02 , + [^^('^1 ,t)-^{r2,t)]-^^ (n , r2 , 0i - 02 , i) 

= -^(ri,r2,0i-02) f-^--/-)c.(ri,tV(r2,t). (C9) 
Introducing the Fourier transforms of the potential of interaction and of the correlation function, we obtain 

— (ri,t) = 2i7r7^ — — / ^2 (ir2 V] n u„(ri, r2)5„(ri, r2, t), (CIO) 
at Ti ovi Jq 



d^ i f \ d 1 (9 \ 

—^{'ri,r2,t) + in p.{ri,t) - rj(r2,t)].g„ = nM„(ri,r2) —- — I u:{ri,t)uj{r2,t). (Cll) 

dt 1 \fi or I r2 or2 ) 

With the Bogoliubov ansatz, Eq. (jCIOp becomes Eq. ([57)) . On the other hand, the asymptotic value of the corre- 
lation function can be obtained by taking the Laplace transform of Eq. ()CIip and considering the limit t — s> +oo. 
Equivalently, it is obtained by making the substitution 

^(ri,r2,t)^ lim e5„(ri,r2,-foo), (CI2) 
dt 6^0+ 

in Eq. (jCTT)) . This yields 

9n{ri,r2,+oo) = — — — — — — — w ri w ra). (C13) 

7 n [ii(ri) — il(r2)\ — ^0+ \ri ori r2 or2 J 

Substituting Eq. (IC13I) in Eq. ()37)) . we obtain Eq. ()45l) which finally leads to the Landau- type equation ([48)) . 
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mentioned as a footnote in his paper); (ii) the inviscid dissipation of energy in 3D turbulence for singular Euler solutions. 

[85] His argument for the existence of negative temperatures was given two years before Purcell & Pound Q] reported the 
presence of negative "spin temperatures" in an experiment on nuclear spin systems. 

[86] An exception concerns the case of non-neutral plasmas under a strong magnetic field 0, U . In the so-called "guiding-center 
approximation" , this system consists in long filaments of charge parallel to the magnetic field B and moving under their 
mutual electric field E with the ExB/B^ velocity. The mathematical description of this system is equivalent to a collection 
of A'^ point vortices in which the charge of a filament corresponds to the circulation of a vortex. 

[87] Actually, the evolution of the distribution function is not due to real collisions between stars (like in the theory of gases) but 
to weak defiexions due to binary encounters with relatively large impact parameter [2^. For simplicity, we shall however 
call them "collisions" . 

[88] Here again, the evolution of the vorticity field is not due to real collisions between point vortices but to distant binary 
encounters. For simplicity, we shall however call them "collisions". 

[89] In fact, it has been proven that the point vortex gas is non-ergodic [sj. However, statistical mechanics does not require 
strict ergodicity. An efficient mixing over the energy surface will suffice to justify the use of the microcanonical ensemble 
and of the Boltzmann distribution for A'^ +oo. However, our concern here is that the point vortex gas may not mix well 
enough due to the effect of "collisions" (correlations, finite A'' effects) to justify the applicability of statistical mechanics 
and of the Boltzmann distribution. 

[90] Contrary to the Landau equation, our derived kinetic equation does not exhibit any divergence, suggesting that collective 
effects are not crucial for point vortices contrary to the case of plasma physics where they account for Debye shielding and 
regularize the large-scale divergence appearing in the Landau equation. 

[91] Nevertheless, a conclusion of our study will be that higher order distributions should also be considered. 

[92] Indeed, we can interprete the kinetic equation (|63p as describing the "collisions" between a test vortex described by the 
variable 1 and field vortices described by the running variable 2. In Eq. (|63p all the vortices are equivalent so that the 
distribution of the field vortices uj{r2,t) changes with time exactly like the distribution of the test vortex uj{r-i,t). In the 
bath approximation, we assume that the distribution of the field vortices uj{r2) is prescribed. 

[93] For point vortices, the kinetic equation (|63|) governing the evolution of the system as a whole does not present any 
divergence contrary to the Lenard-Balescu equation in plasma physics that presents a divergence at small scales. However, 
a logarithmic divergence at small scales occurs in the Fokker-Planck equation (|66|) when we make the bath approximation 
and consider that the dominant contribution to the (5-function comes from the interaction at r' = r. This divergence comes 
from the assumption made in the kinetic theory that the point vortices essentially follow the streamlines produced by the 
mean flow (modifled by collective effects). This neglects the contribution of "hard" collisions with small impact parameter 
that lead to more complex trajectories and more complex interactions. As explained in Sec. Ill Bl these hard collisions 
could be taken into account by keeping the contribution of the fourth term in Eq. (|17|l . Alternatively, we can proceed 
heuristicaUy and reg ularize the logarithmic divergence by introducing a cut-off at the scale at which hard collisions come 
into play [39l. [SOl. l6l| . With this regularization, it can be shown that InA ~ ^InA'^, leading to a InAf factor in the Fokker- 
Planck equation (|66|l . Since only distant collisions r2 7^ ri occur in the kinetic equation (|63|) describing the evolution of 
the system as a whole, there is no divergence at close separation, and therefore no In A'' factor, in that equation. 



